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A general formulation of the Boundary Quantum Inverse Scattering Method 
is given which is applicable in cases where i?-matrix solutions of the Yang- 
Baxter equation do not have the property of crossing unitarity. Suitably modi- 
4h . fied forms of the reflection equations are presented which permit the construc- 

tion of a family of commuting transfer matrices. As an example, we apply the 
i -^h ■ formalism to determine the most general solutions of the reflection equations for 

a solution of the Yang-Baxter equation with underlying symmetry given by the 
Drinfeld double D{D^) of the dihedral group D%. This i?-matrix does not have 
the crossing unitarity property. In this manner we derive integrable boundary 
conditions for an open chain model of interacting non-abelian anyons. 

> 
OO 

cti ' 1 Introduction 

The study of systems with non-abelian anyonic degrees of freedom currently attracts 
high interest, due to the possibilities for exploiting their topological properties to 
encode quantum information in a manner which is protected from decoherence pQ. 
j> | An appropriate framework in which to formulate systems with anyonic symmetries 

is through the representation theory of quasi-triangular Hopf algberas [2], which in- 
cludes the class of Drinfeld doubles of finite group algebras [3j H] . This latter class 
of algebras is particularly suited for the description of non-abelian anyons where the 
conjugacy classes and centraliser subgroups of the finite group label generalised no- 
tions of the magnetic and electric charges [HE]- Moreover, within the quasi-triangular 
Hopf algebra framework, consistent braiding and fusion properties for anyonic the- 
ories are naturally obtained. The braiding properties are characterised by solutions 
of the Yang-Baxter equation without spectral parameter, which are realised through 
the universal .R-matrix of the algebra. The fusion properties are given by decompo- 
sitions of tensor product representations of the Hopf algebra, which are governed by 
the coproduct structure. These fusion rules provide a means to construct interacting 
systems by assigning energies to the various possible multiplet structures. This then 
enables the study of one-dimensional (chain) models with local interactions, as has 
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been recently undertaken in [7j using Fibonacci anyons. For this case the interaction 
energies were chosen in such a way that the local Hamiltonians provided represen- 
tations of the Temperley-Lieb algebra, which necessarily means that the system is 
integrable and can be solved exactly. A study of a non-integrable non-abelian anyon 
chain can be found in [Sj. 

The theory of integrable chains has a long history associated with the Quantum 
Inverse Scattering Method (QISM) [9], which relies on a solution of the Yang-Baxter 
equation with spectral parameter to construct a family of commuting transfer ma- 
trices. The transfer matrix may be used to generate the conserved operators of an 
integrable quantum system. Following this procedure we have previously shown that, 
using the Drinfeld double D(D 3 ) of the dihedral group P 3 , there exists a spectral 
parameter dependent solution of the Yang-Baxter equation which can be used to 
construct an integrable interacting non-abelian anyon chain 1 1 j . There the standard 
approach of the QISM was used, producing a chain with periodic boundary condi- 
tions. For open chain cases integrable boundary conditions are provided by solutions 
of the reflection equations, as was first elucidated by Sklyanin [TTJ, and is generally 
known as the Boundary Quantum Inverse Scattering Method (BQISM). Our goal is 
to extend the BQISM formalism in a manner which will enable the construction of a 
non-abelian anyon open chain with integrable boundary conditions. 

In Sklyanin's original formulation of the BQISM several conditions were imposed 
on the P-matrix including P-symmetry, T-symmetry and crossing symmetry [TTJ. 
It was soon realised that the BQISM can be extended to cases where the P- and 
T-symmetry properties are relaxed to the more general PT-symmetry |12j . and the 
crossing symmetry property can be replaced by the more general crossing unitarity 
condition [T3] (see equation ffTUj) below for the definition). Later it was shown in [T3] 
that the BQISM can be formulated for cases without PT-symmetry. Here we further 
extend the formulation of the BQISM by removing the imposition of the crossing 
unitarity property. This is necessary to construct integrable boundary conditions for 
the D(D$) anyon chain, as the P- matrix does not possess this property. 

In Section 2 we present the formulation of the BQISM for P-matrices without 
crossing unitarity. Using the explicit example provided by the D(D 3 ) P-matrix of |lUj . 
in Section 3 we explicitly find the most general solutions of the reflection equations. 
In Section 4 we use these results to derive a non-abelian anyon chain with integrable 
boundary conditions, and concluding remarks are given in Section 5. 

2 BQISM for i?-matrices without crossing unitar- 
ity 

Our first objective is to reformulate the BQISM with a minimum number of assumed 
properties imposed on the P-matrix solution of the Yang-Baxter equation. We start 
with invertible operators R(z) G End {V ®V) and L(z) G End {V ®W) which satisfy 
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the Yang-Baxter equation on End (V <8> V <8> V): 

Ri2(xy- l )R n (x)R 2Z (y) = R 23 (y)Ri 3 (x)R 12 (xy- 1 ), (1) 
and the intertwining relation on End (V <S> V" <S> W) : 

R 12 (xy- 1 )L 13 (x)L 23 (y) = L 23 (y)L 13 (x)R 12 (xy- 1 ). (2) 

Here the subscripts indicate on which vector spaces each operator acts, so for example 
^23(2) = I ® -R(z). Solutions to the Yang-Baxter equation (fj]) are referred to as R- 
matrices, while L(x) appearing in (j2J) is called an L-operator. We impose only the 
following conditions on R(z): 

1. R tl {z) is invertible and 

2. obeys regularity, i.e. R(l) = P- 

Here t\ denotes the partial transpose over the first space and P G End (V ®V) is the 
usual permutation operator defined by 

P(v ® w) — w ® v , v, w e V". 

The following theorem is reproduced from |15j : 

Theorem 2.1. If R(z) is an R-matrix satisfying the regularity property then it also 
satisfies the unitarity property, i. e. 

R l2 {z)R 2l {z- 1 ) = f(z)I®I, 

with f(z) a scalar function satisfying f(z) = f(z~ l ). 

Proof. Let R(z) be an .R-matrix satisfying regularity. Then 

R 12 (z)R 13 (l)R 23 (z- 1 ) = R^z-^R^Ruiz) 

=> R l2 {z)P l3 R 23 {z- 1 ) = R 23 {z' l )P l3 R l2 {z) 

=> i? 12 (^) J R 21 (^ 1 )P 13 = R 2 z{z- l )R 32 (z)P l3 

R l2 {z)R 2l {z- v ) = R 23 {z~ v )R 32 {z) 

The left hand side acts trivially on the the third space while the right hand side acts 
trivially on the first. Combining these it follows that R\ 2 {z)R 2 i{z~ l ) must be a scalar 
of the identity, and that the scalar function is invariant under z — > z . □ 

Utilising the condition that R tx (z) is invertible we define the operator 

n 12 (z) = [(RlMzTT, (3) 

which by definition implies 

nUz)R^(z) = n x {z)R%{z) = 1. (4) 
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We now introduce two reflection equations 

R 12 (xy- 1 )K^(x)R 21 (xy)K^(y) = K 2 ~ (y)R 12 (xy)K^ (x)R 21 (xy- 1 ) J (5) 
R 12 (yx- 1 )K+(x)TZ 21 (xy)K+(y) = Ki{y)TZ 12 {xy)Kt{x)R 21 {yx- 1 ), (6) 

where K + (z), K~(z) G End (V) are known as reflection matrices. The matrices 
K + (z), K~(z), R(z) and L(z) will enable us to construct an integrable model on an 
open chain. The transfer matrix is defined as 

t(z)=tr a [K+{z)T(z)], 

where tr a is the trace over space a and T(z) is the double monodromy matrix 

T(z) = L aN (z)...L al (z)K-(z)L- 1 1 (z- 1 )...L^(z- 1 ). 

It is known and easily verifiable from §2§ and (j3J) that 

Ri2(xy~ 1 )Ti 3 (x)R 21 (xy)T 23 (y) = T 23 (y)R 12 (xy)Ti 3 (x)R 21 (xy~ 1 ). (7) 

Proposition 2.2. The transfer matrices t(x),t(y) commute for all x,y e C. 

Proof. 

fixy-'XxXy) 

= (tn ®tr 2 ) {fixy-^K+MK+ix^TilWTM} 

= (tn ® tr 2 ) {/(xy- 1 )ir 2 + (y)ir i + (a;) t ^i 1 2 (xy) J R 2 1 1 (xy)T 1 ^(x)T 2 3(y)} 

= (trx ® tr 2 ) {/(^- 1 )ir+(y)^ 12 (^)ir i + (x)T 13 (a;) J R 21 (a;y)T 2 3(y)} 

= (tn (g) tr 2 ) { [ J fr+(y)7e 12 (^)ii: i + (x)i2 21 (yx- 1 )] [R^xy-^T^R^xy)^)] } 

= (tn ® tr 2 ) { [ifofox - Wfr^iM^l/)] [T 23 (y)i? 12 (a; 2 /)T 1 3(x) J R 21 (xy- 1 )] } 

= (tri ® tr 2 ) {/(xy- 1 )ir i + (x)^ 21 (xy)^+(y)T 2 3(y) J R 12 (x 2 /)T 13 (x)} 

= (tn ®tr 2 ) {/(xy- 1 )^^)^^)* 2 ^^)^^)^^)^^)} 

= (tn ® tr 2 ) {/(^- 1 )K 1 + (x)K 2 + ( 2 /)* 2 T 2 ^( 2/ )T 13 (x)} 

= f{xy- l )t{y)t{x) 

where we have used equations PfoTH) and Theorem 12.11 □ 
We now impose the limit condition 

K-(1)=I 

and only consider instances for which tr(i^ + (l)) ^ 0. Then in the case L(z) = 
R(z) the global Hamiltonian on an open chain with boundary fields is defined in the 
following way: 



H 



2 = 1 



2tr(AT+(l)) [dz 

N-l , 

1=1 



+ 



2=1 



Z=1J 

tTg (C^ljgjvj 

tr(tf+(l)) ' 



where c G C and the local Hamiltonians are given by 



2=1 



We will refer to the first term of ([H]) as the bulk Hamiltonian, while the second and 
third terms describe boundary field interactions. By construction the global Hamilto- 
nian (E]) commutes with the transfer matrix t(z), which means that the Hamiltonian 
is necessarily integrable. The conserved operators commuting with the Hamiltonian 
are obtained as the co-efficient operators in the series expansion of t(z). 



3 Reflection matrices for an i?-matrix associated 
with D(D 3 ) 

We now apply the above formalism to solve for the reflection matrices satisfying 
equations (0) and This R- matrix we use is constructed from the representation 
theory of D(D^) [10], with associated L-operators given in [16]. Explicitly, we have 
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(9) 



The properties of the .R-matrix include regularity, and consequently unitarity. How- 
ever it can be verified that the i?-matrix does not satisfy the crossing unitarity con- 
dition, i.e. there does not exist M G End(V) and A G C such that 



(10) 



This is in contrast to i?-matrices obtained from loop representations of affine quantum 
algebras, for which equation (TTTJT) is always satisfied [T7] . 
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We first calculate H(z) defined by equation (jHJ): 
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These are the two operators required to construct the reflection matrices K (z) 
K + (z). 



and 



3.1 Special case of the reflection equation 

To determine the possible matrices, K~(z) and K + (z), which satisfy equations §5§ 
and (El) it is first convenient to determine all non-diagonal invertible matrices K(z) 
which satisfy the equation 

K 2 (y)R 12 (0)K 2 (z Q )R 12 (0) = R 12 (0)K 2 (z )R l2 (0)K 2 (y), (11) 

where y,z G C, z is fixed and -R(O) = PR(0). We scale K(z) so that the entries 
of lim 2 ^ 2o K(z) are all finite and that at least one is non-zero, as is always possible. 
Throughout this section we write K(z) in the form 

3 

K{z) = KAz)E) 

where Ej denotes the elementary matrix with a 1 in the ith row and jth column. We 
consider the indices of the functions h%j{z) and elementary matrices Ej modulo 3. 
Using this notation, K(z) is a solution to equation (TTTj) if and only if 

h id (z )h k j(y) = hij +k+ 2i(z )h 2 i + 2k,2i+2i(y) 
for all 1 < i, j, k, I < 3. 

Proposition 3.1. If K(z) satisfies equation (1 1 il) and h a ^ a (zo) = for some 1 < a < 3 
then h a j(z ) = for all 1 < j < 3. 

Proof. Assume there is an integer b such that h a ^(z ) ^ 0. As K(z) satisfies equation 
(TTTj) . we have 

h a ,a( Z o)ha,l(y) = {zq) h a ,2a+2l (y) h ajb (z )ha,b(y) = 0. 

This contradiction proves the proposition. □ 
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Corollary 3.2. If K(z) is a solution to equation f TTTT) then K(zq) has at least one 
non-zero diagonal entry. 

Proof. Assume there is a solution where all the diagonal entries of K(zq) are zero. It 
follows from Proposition 13.11 that h at j(z ) = for all 1 < a,j < 3, which contradicts 
our requirement that at least one entry of K(zq) is non-zero. □ 

Proposition 3.3. If K(z) is a non- diagonal matrix satisfying equation 07]) and 
h a , a ( z o) 7^ for some 1 < a < 3 then h aa+b (z ) ^ for b G {1,2}. Furthermore, 
h a ,a+2{z) = ah a>a+1 (z) and a 3 = 1 where a = h a , a+1 (z ) / h a , a (z ) . 

Proof. We now assume that there exists a b G {1,2} such that h ata+ b(z ) = 0. This 
leads to 

h a ,a+b(zo)K,a+b(y) = K,a( Z o)h a ,a+2b(y) => K,a+2b(y) = 0. 

Therefore h a ^ a+2 b(y) = 0, which implies through the same argument that h a a+b (y) = 
0. Hence 

h a ,a+l(y) = h a , a +2(y) = 0. 

But 

h a ,a(Zo)hl +c ,l(y) = h^ a+c (z )h 2a+2 l + 2c,2a+2l{y) => hi+c,l{y) = 

for 1 < / < 3 and c G {1,2}. This implies that if h ata+ b(z ) = for some integer 
b G {1,2} then K(z) is diagonal, which is a contradiction. Hence h aa+b (z ) ^ for 
fee {1,2}. 

To show the other half of the proposition we set a = h a ^ a+1 (z ) / h a a (z ) . Then 

(y) = h a>a (zo)h a>a+2 (y) =>- h a>a+2 (y) = ah a)a+1 (y) and 
ha,, a +2{za)h a>a+2 (y) = h a , a (z )h a , a+1 (y) [a 3 - l]h a , a +i(y) = 0. 

This completes the proof. □ 

Proposition 3.4. For a non-diagonal matrix, K(z), which satisfies equation 177]) 
and has h a ^ a (zo) ^ and h a+ i^ a+ i(z ) ^ for some 1 < a < 3 then h a+ 2 )a +2(zo) ^ 0. 

Proof. We find 

K+l,a+l(zo)h a+ 2,a+2(y) = h a+lja+1 (z )h ata (y) h a+2 ,a+2(y) = K,a(y) i 

which proves the proposition. □ 

Proposition 3.5. Any non-diagonal matrix, K(z), which satisfies equation ( 1771) and 
has hi^(zo) 7^ for all 1 < i < 3 is of the form 

/ A(z) aB(z) a 2 B(z) \ 
K(z) = P 2 B(z) A(z) (3B(z) 

V jB(z) j»B(z) A(z) J 

where a 3 = (3 3 = 7 3 = a(3 r ) = 1 and A(zq) = B(zq) = 1. 
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Proof. As hx t i(zo) ^ and we are free to scale by a constant we set hi t i(zo) = 1. We 
find that 

hi,i{z )[hi t i(y) - h 2i+ 2i,2i+2i(y)] =0 =>> = 
for all 1 < i, j, I < 3. 

We now define the following variables 

A a = h a ^ a+ i(zo), 

for 1 < a < 3. We see that from the previous proposition that A^ = 1 and 

K,a+2iy) = KK,a+l(y) 

for all a G {1,2,3}. This gives the properties required for within each row, but we 
still need to relate the entries down each column. We have that 

h a ,a+c{Zo)h a+ 2c,a+ '< / a,a+2c 

for a G {1, 2, 3} and c G {1, 2}. Hence 

(y) = KK+i, a +2(y) =>- A a+2 /i a+2>a (?/) = A a /i a+ i ia+2 (y) 

Expressing this explicitly we have that 

AiA 3 /ii, 2 (y) = A 2 Ai/i 2 ,3(y) = A 2 A 3 /i 3i i(y) 

Hence the off-diagonal entries are all scalar multiples of each other. Moreover, by 
considering the above equation at y = zq we find 

A^A 3 = AiA 2 AiA 3 A 2 = 1. 

This proves the proposition. □ 

Proposition 3.6. If K(z) satisfies equation H7\) and K(zq) has only one non-zero 
diagonal entry then K(z) can be written in the form (after basis transformation and 
scaling) 

/ A{z) aB{z) a 2 B(z) \ 
K(z) = aD(z) C(z) E{z) 

\ D{z) aE{z) C(z) J 

where a 3 = 1, A(z Q ) = B(z ) = 1 and C(z ) = D(z ) = E(z ) = 0. 

Proof. First note that R is invariant under relabelling of the indices, so without loss 
of generality we can set hi t i(zo) = 1- We set A(z) = hi ; x(z), so A(z ) = 1. By 
Proposition 13.31 we know that /ii i2 (^o) ^ 0, and thus we set /ti, 2 (z) = ctB{z) with 
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B(zq) = 1. By Proposition 13.31 we also have that a 3 = 1 and ^1,3(2) = a 2 B(z). 



We now consider the diagonal entries and see that 

hi,i{z )h 3t3 (y) = /n,i {z )h 2>2 (y) ■ 
It follows that h^ 3 (z) = h 2 ^ 2 (z). We let C(z) = h 2 , 2 {z). 
We now use 

hi,i(z )h 2 ,i(y) = h 1)2 i(za)h 3i 2+2i{y)- 

This gives 

h 2 ,i(y) = ah 3A (y) and h 2 , 3 (y) = a 2 h 3i2 (y). 

We let D(z) = h 3< i(z) and E(z) = h 2j3 (z). Lastly, note that C(z ) = and hence by 
Proposition 13.11 we have D(z ) = E(z ) = 0. □ 

We have classified all the possible non-diagonal matrix solutions to equation f fTTT) . 
These can be classified by the number of non-zero diagonal elements of K(zq). In 
Corollary 13.21 and Proposition 13.41 it was shown that K(z) can only have one or three 
non-zero diagonal elements. Moreover, these solutions can all be written (after basis 
transformation and scaling) in the forms given in Propositions l3.5l and l3.6l respectively. 

3.2 Reflection matrix K~(z) 

We now turn our attention to the reflection equation ©. Setting x = 00 and using 
regularity, we obtain the equation 

Ra(0)K^ (00)^x2(0)^- (y) = K 2 ~(y)R 21 (0)K{(oo)R 12 (0). 

This is equivalent to equation f fTTT) with zq = 00. Using the previous section we know 
all the possible forms of K~(z). We now only need check under which conditions, if 
any, these forms satisfy the reflection equation (J5j). 

Proposition 3.7. The only diagonal matrices satisfying equation |3]] are scalars of 
the identity. 

Proof. The proof involves a straightforward calculation so we omit the details. □ 
Proposition 3.8. There are no invertible matrices satisfying equation |3P of the form 

f A(z) aB(z) a 2 B(z) \ 
K'(z) = aD{z) C(z) E(z) 

\ D{z) aE{z) C(z) J 

where a 3 = 1, A(oo) = B(oo) = 1 and C(oo) = D(oo) = E(oo) = 0. 
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Proof. We first assume that such a K (z) does exist. Substituting K (z) into the 
reflection equation, we find it must satisfy the constraints 

D(y)B(x) = B(y)D(x) and E(y)B(x) = B(y)E(x). 

This implies that D(z) and E(z) are both scalars of B(z). Using the boundary 
conditions at z = oo we deduce that 

D(z) = E(z) = 0. 

After imposing D(z) = E(z) = we obtain the second constraint equation 

B(x)B(y) = 0. 

This is a contradiction, as B(oo) = 1, and hence there are no solutions K~(z) of the 
above form. □ 

Proposition 3.9. All matrices satisfying equation ([3]) are either scalars of the iden- 
tity or of the form 

w 2 + bz- z 2 (1 - z 2 ) (1 - z 2 ) 

K - ( ) = | w 2 (1 - z 2 ) w 2 + bz-z 2 w(l-z 2 ) 
w 2 (l — z 2 ) w(l — z 2 ) w 2 + bz — z 2 

for some b G C and w a cube root of unity, up to scaling and a change of basis. 

Proof. Let K~{z) be a solution to equation ([5]) which is not a scalar of the identity. 
Then by the earlier propositions we know K~(z) can be written in the form 

A(z) aB(z) a 2 B(z) 
K~(z) = ( (3 2 B(z) A(z) (3B(z) 
1 B(z) i 2 B(z) A(z) 

where a 3 = j3 3 = 7 3 = a(3^ = 1 and A(oo) = B(oo) = 1. Substituting this into 
equation (0), we obtain the constraint equation 

= -x(x 2 + a 2 f3x 2 + a/3 2 ) - x(x 2 + aft 2 + a 2 (3)y 2 
+ (a[3 2 x 2 + (3 2 a + x 2 + x 4 + 2a 2 f3x 2 )y 
-3/(1 - x 2 ) 2 f(x) + x(x 2 - l)(y 2 - l)f(y) 

where f(z) = A(z)/B(z). As the coefficient of f(x) is linear in y, we take the double 
derivative with respect to y, obtaining 

o (x 2 + a(3 2 + a 2 (3) , 2 

2^ -r- 2 tt = 2/(3/) + 4yf'(y) + (y 2 - 1 / {y . 
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The right hand side is independent of x, so we deduce that 

= 1 + a(3 2 + a 2 (3 = 1 + a 2 (3 + (a 2 p) 2 . 
This implies that a 2 /3 is a primitive cube root of unity. We set 

(3 = wa, 

where w is a primitive cube root of unity. Our matrix becomes 

/ A(z) aB(z) a 2 B(z) 
K~(z) = w 2 a 2 B(z) A(z) waB(z) 
\ w 2 aB(z) wa 2 B(z) A(z) 

The different choices of a are equivalent up to a basis transformation under which 
R(z) is invariant, so without loss of generality we choose a = 1. Now the differential 
equation reduces to 

2 = 2f(y)+4yf(y) + [y 2 -l}f"(y). 
This has the general solution 

f( n a + by-y 2 
f{y) = — ; — 

1 - y 2 

where a, b G C. Without loss of generality we set 

1 - z 2 

A(z) = 1 and B(z) 



a + bz — z 2 



Substituting K (z) into equation (jHJ), we find that we must have a = w 2 . Furthermore 
if a = w 2 then K~(z) satisfies equation (jSJ) for all b G C. □ 

3.3 Reflection matrix K + (z) 

In this section we construct solutions to the other reflection equation (jBJ). Setting 
x = and using unitarity, we obtain the equation 

R 21 (0)K+(0)K 21 (0)K+(y) = K+(y)K 12 (0)K+(0)R 12 (0), 

This is equivalent to equation (iTTj) with ^ = as 7^ 2 i(0) = i?i2(0). Using the previous 
section we know all the possible forms of K + (z). It remains to check whether these 
forms provide solutions to the reflection equation and if so, under what conditions. 

Proposition 3.10. The only diagonal reflection matrices satisfying equation (0|) are 
scalars of the identity. 

Proof. The proof involves a straightforward calculation so we omit the details. □ 



11 



Proposition 3.11. There are no invertible matrices satisfying equation (0|) of the 
form 

/ A(z) aB(z) a 2 B(z) \ 
K + (z) = aD(z) C(z) E(z) 

\ D{z) aE{z) C{z) j 

where a 3 = I, A(0) = B(0) = 1 and C(0) = D(0) = E(0) = 0. 

Proof. Directly substituting K + (z) into equation we find the constraints 

D(y)B(x) = B(y)D(x). 
This implies that D(z) is a scalar of B(z). As D(0) =0 we deduce that 

D{z) = 0. 

We also obtain the constraint equations. 

C(y)A(x) = A(y)C(x) and E(y)B(x) = B(y)E(x) 
Using similar reasoning, we find that 

E{z) = C{z) = 0. 

Hence K + (z) is not invertible, which proves the proposition. □ 

Proposition 3.12. All matrices satisfying equation are either scalars of the iden- 
tity or scalars of 

(1 + bz- wz 2 (1 - w 2 z 2 ) (1 - w 2 z 2 ) \ 

w 2 (l-w 2 z 2 ) 1 + bz-wz 2 w(l-w 2 z 2 ) , 

w 2 (l — w 2 z 2 ) w(l — w 2 z 2 ) 1 + bz — wz 2 J 

for some b G C and w a primitive cube root of unity. 

Proof. By the earlier propositions, if K + (z) is not a scalar of the identity then it must 
be of the form 

/ A(z) aB(z) a 2 B(z) \ 
K + {z) = P 2 B(z) A(z) (3B(z) 

V iB(z) 7 2 B(z) A{z) ) 

where a 3 = /3 3 = 7 3 = a(3^ = 1 and A(0) = B(0) = 1. Using a similar technique 
to that used to prove Proposition 13.91 we find a constraint equation and differentiate 
twice with respect to x. We obtain the following two equations: 

2a/? 2 = 2f(y) + Ayf'(y) + [y 2 + l + a(3 2 ]f"(y) 
= (l + pa 2 + (Pa 2 ) 2 )[2f'(y)+yf"(y)] 

where f(y) = A(y)/B(y). These equations only have a solution in common when 

(1 + (3a 2 + ((3a 2 ) 2 ) = 0. 
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Hence we must have (3 = wa = w 2 ^ where w is a primitive cube root of unity. Fur- 
thermore, the three different choices of a are equivalent up to a basis transformation 
under which R(z) and lZ(z) are both invariant, so without loss of generality we can 
choose a — 1. The differential equation simplifies to 

2w 2 = 2f(y) + Ayf (y) + [y 2 - w]f"(y). 

This has the general solution 

a + bz — wz 2 

f{z) = — 

1 — w z z z 

where a, b G C. However /(0) = A(0)/S(0) = 1, so a = 1. We set 

= 1 + bz - wz 2 and B(z) = \-w 2 z 2 . 

Substituting this into the reflection equation ([6]), we find this is a solution for all 
b G C. This completes the proof. □ 



An integrable Hamiltonian with open boundary 
conditions 



We found in the previous section all possible solutions K + (z) and K~(z) of the 
reflection equations (jSJ) and ([6]) for our i?-matrix with non-abelian anyonic symmetry. 
To calculate the global Hamiltonian with non-trivial boundary terms, we use the 
reflection matrices in the following form: 



/ 



K~{z) 



1 



aw 



w^+az—z 2 
w 2 +az — z 2 



a 2 w 



w 2 +az— z 2 
1 

(I-* 2 ) 
ui 2 +az— z 2 



(l-^ 2 ) \ 

.9 , , -.9 



aw 



«> 2 +az— z 2 
(l-^ 2 ) 



«; 2 +az- 
1 



and 



JT+(s) 



1 + bz — w j z 2 
w 2 ip 2 (l 
w 2 if3{l 



w 2 h 2 ) 
w 2 iz 2 ) 



(3(1- 
1 + bz 
wi(3 2 (l 



w 2j z 2 



) /? 2 (1 



w 



Z 2 ) 

w j (3(l - w 2j z 2 ) 
w 2j z 2 ) 1 + bz — w j z 2 



w j z 2 



where a, ft G {1, to, w 2 }, j G {1, 2}, a, 6 G C and w is a primitive cube root of unity. 
As we are only considering the case where tr(i^ + (l)) ^Owe impose that b ^ w j — 1. 
Choosing c = i in (jSJ) yields the local Hamiltonian as given in [H 



P7(l) 
V 7(2) 



7(2) 
7(3) 



E 



7(2) 
7(3) 



E 



, 7 (1) N 
7(2)' 



£ «<3 

7e-D 3 

where the elements 7 G -D3 are written as permutations of {1,2,3}. (Recall that 
D3 is isomorphic to the permutation group S3.) The bulk Hamiltonian, which de- 
scribes nearest neighbour interactions between anyons with local three-dimensional 
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state spaces, commutes with the action of D(D 3 ) algebra [TP) . For the boundary 
terms we find 

5 1™ = A 



and 



tr 8 (iJif+(l)^, a ) =B 






a 


a 2 \ 




n 

u 




aw 2 


a 2 w 


o / 





/? 




2 ^' 







to; 3 ' 









-w j + b)- 1 e 



tr(Jf+(l)) 

where A = z(l — w; 2 + a)~ x G C and Z? = — w 3 + b)^ 1 G C. We choose a 
iuX^ 1 + cj -1 — 1, b = w J — 1 — iu)~ 3 Y~ x with X, Y G K unconstrained parameters. 
This leads to A = B = w 3 Y, resulting in a global Hamiltonian which is 

hermitian. 

Alternatively, one can choose the reflection matrices to be K + (z) = K~(z) = I. 
This yields an open chain Hamiltonian without boundary interaction terms. It is also 
possible to have one end of the chain with a boundary interaction term, while the 
other end is without a boundary interaction term. 



5 Summary 

We reformulated the BQISM through a pair of reflection equations in a fashion which 
does not rely on the i?-matrix solution of the Yang-Baxter equation satisfying the 
crossing unitarity condition. This was motivated by the case of the i?-matrix as- 
sociated with the quasi-triangular Hopf algebra D(Ds), for which crossing unitarity 
does not hold. We then proceeded to determine the most general solutions of the re- 
flection equations. With these results we were able to determine integrable boundary 
conditions for an anyonic chain where the bulk Hamiltonian has D(D 3 ) symmetry. 

The i?-matrix associated with D(D^) is the simplest example in a heirarchy of 
solutions associated with D(D n ) which solve the Yang-Baxter equation [T8J. All of 
these solutions are characterised by an absence of crossing unitarity, meaning that 
our general formalism is applicable on a wider scale. 

A future direction in this program of research is to compute the Bethe ansatz so- 
lution of the open chain model derived above. Implementation of the algebraic Bethe 
ansatz appears problematic in this case, due to the lack of a suitable pseudovacuum 
state. A more promising avenue is offered by a functional approach aided by fusion 
relations e.g. see [Tj5] and references therein. 
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